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Abstract. Low-temperature screening properties of the inhomogeneous two-dimensional 
electron gas in a Hall bar subjected to a strong perpendicular magnetic field are explored using 
a self-consistent approach. An external oscillating modulation potential with an amplitude of 
the order of the cyclotron energy is added to the electron-confining background potential, and 
the resulting change of the self-consistent potential is investigated as function of modulation 
strength, magnetic field, and temperature. The consequences of Landau-level pinning and 
the interplay of compressible and incompressible regions for the resulting strongly non-linear 
screening phenomena are explained. 



1. Introduction 

As a consequence of the highly degenerate, Landau-quantised energy levels, a two- 
dimensional electron gas (2DEG) in a strong perpendicular magnetic field has unusual low- 
temperature screening properties |[T1I2]|. In an inhomogeneous 2DEG with long-range density 
fluctuations, pinning of Landau levels (LLs) leads to quasi metallic regions with high density 
of states (DOS) at the Fermi energy and therefore nearly perfect screening ability. These 
so called "compressible" regions coexist with quasi-insulating "incompressible" regions 
separating compressible regions with different LLs at the Fermi energy. In the incompressible 
regions the Fermi energy falls into the gap between two adjacent LLs and, in thermodynamic 
equilibrium at low temperatures, no redistribution of electrons is possible which could 
contribute to screening. Also, in the incompressible regions the electron density nei{r) is 
constant, corresponding to total filling of an integer number of LLs, while in the compressible 
regions nei(r) adjusts itself so that the self-consistent electrostatic potential energy V{r) of 
an electron differs from the Fermi energy, more precisely from the constant electrochemical 
potential /i* by a Landau energy fkJc^n + 1/2), where cUc = eB/m is the cyclotron frequency 
in the magnetic field B. As a consequence, V{r) differs between different compressible 
regions by integer multiples of HiOc- Landau level pinning and the interplay of compressible 
and incompressible regions will strongly affect the screening properties of the 2DEG. 

In a bounded 2D geometry with translation invariance in one direction, the compressible 
and incompressible regions degenerate to strips parallel to the boundary. For half-space and 
Hall-bar geometries models with planar charge distributions have been proposed that allow 
closed solutions of Poisson's equation, i.e., the calculation of the potential for given electron 
density, and estimates of position and widths of the incompressible strips have been given 
0|. By adding the non-linear Thomas-Fermi approximation for the calculation of the 
electron density from the potential, that work was extended to a self-consistent approach 
which allows to calculate both electron density and electrostatic potential for arbitrary 
temperature [5, 6|. This approach shows that the existence and the width of incompressible 
strips depends sensitively on temperature. It also permits to calculate their position and width 
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for given background charges without additional assumptions. We adopt this approach to 
investigate systematically the screening of a harmonic external potential in such a confined 
2DEG. 

2. The Thomas-Fermi-Poisson (TFP) approach 

In order to employ the approach presented in Ref. f6l, we assume that our 2DEG lies in the 
2; = plane, laterally confined by in-plane gates located dXx < —d and x > d, which are kept 
at voltages Vl and Vr, respectively. In this work we assume a symmetric depletion on both 
sides with the depleted strips having a width of \d — b\, and we will consider only symmetric 
potentials V{—x) = V{x), notably with V{±d) = Vl = Vr = 0, which leave the electron 
density profile symmetric, nei{—x) = nei{x). Then, solution of Poisson's equation yields for 
the electrostatic potential energy of an electron 



(1) 



2e r 

V{x) = — dtK{x,t) p{t), K{x,t)=\n 

K [X — t)d 

where p{x) is the surface charge density in the strip \x\ < d and k is an average dielectric 
constant 1 6 1 . If we assume a homogeneous surface density hq of positive background charges, 

p{x) = e[no — nei{x)], we obtain 

V{x) = V,g{x) + VH{x), where V,g{x) = -E^^l - {x/ df (2) 

with Eq = 27re^nQd/Ris the confining potential due to the background and Vf/(x) the Hartree 
contribution due to the electron density, given by Eq. ([T]) with p(x) replaced by —enei{x). 
Assuming that V{x) varies on a characteristic length much larger than typical quantum 
lengths, notably the magnetic length = a/ (h/eB), one can calculate nei{x) for given 
V{x) in the Thomas-Fermi approximation |l6|| 

n,i{x) = j dED{E)f{[E + V{x) - /i^/fc^T), with /(e) = [1 + e^' (3) 

the Fermi function, p* the electrochemical potential and T the temperature. Here D{E) is the 
density of states of the 2DEG, which we take to be the bare spin-degenerate Landau DOS, 

00 

D{E) = Yl - ^^(^ + 1/2)). (4) 

n=Q 

This completes the TFP scheme for the self-consistent calculation of n(,i{x) and V{x) for 
given temperature and magnetic field and homogeneous background charge. To study 
screening effects, we add to Vhg{x) in Eq. ^ the harmonic external potential 

Vm{x) = Vocos{kxx), with kx = {\ + l/2)Tx/d, (5) 

where A is an integer (to preserve the boundary conditions), calculate the resulting self- 
consistent potential and electron density as function of the modulation strength Vq, and 
compare these results with those for Vq = 0. [In principle Vm{x) can be generated by a 
redistribution of the background charge density.] 

For technical reasons we always start the iterative solution of the self-consistent TFP 
scheme with a calculation for vanishing temperature and magnetic field and for homogeneous 
background, i.e., we take nei{x) = Do[Ef — V{x)]6{Ef — V{x)) with Dq = m/nh? as 2D 
DOS and define the symmetric density profile by the requirement V{b) = V{—b) = Ep = 
p*{T = 0). Then Eq. dU) reduces to a linear integral equation. 



V{x) = -^0^1-^)' + ^ //^f^^ ~ ^^^^^ ^^"^'^^ 
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with oo = fi:/27re^L'o the bare screening length. [Due to a misprint, in the corresponding 
Eq. (15) of Ref. the factor vr in the denominator is missing.] The choice of h defines the 
average density, which is then kept fixed for the calculations at finite temperature, magnetic 
field, and modulation. 

3. Results and discussion 

In the numerical calculations we keep some parameters fixed. So we divide the interval 
—d < X < d into 500 subintervals and calculate V{x) and nei{x) on the corresponding 
equidistant mesh points x = x„. The depletion length is defined by taking b/d = 0.9 and the 
bar width by nao/d = 0.01. For GaAs values (R = 12.4) uq ~ 5nm, this means 2d ^ 3 /xm, 
which is small compared to typical Hall bars used in experiments. Calculation for smaller 
a^/d, i.e., larger d requires more subintervals, and thus larger storage and computation time, 
to achieve the same accuracy near the incompressible strips, but yields no qualitative changes. 
Therefore, we fix the mentioned parameter values and obtain from the solution of the system 
of linear equations defined by Eq. Q an average electron density that yields the average LL 
filling factor z/ = 2 for a cyclotron energy fkUc = = 0.2311 x 10~^-E'o- Starting from this 
solution we find the solution of the non-linear TFP scheme at finite T and B implementing the 
Newton-Raphson method Ollll. For the bare modulation potential © we take A = 2, so that 
Vm{x) has a maximum at a; = and two minima at x = ±OAd, i.e. well inside the Hall bar, 
whereas near the maxima at x = ±0.8d the confinement potential increases strongly and the 
electron density starts to decrease, so that the screening properties are obscured by boundary 
effects. 

Figure[T]shows the potentials and electron densities of the Hall bar calculated for several 
modulation amplitudes at strong magnetic fields and low temperatures. On the left panel 
we consider the case Vq = and show the temperature dependence in the inset, while 
the right panel shows the evolution of the potential and electron density with increasing 
modulation amplitude. All the electron densities are expressed in terms of local filling given 
by z/(s) = 2Txn(,i{x)l'^ and potentials are normalised by Eq. 

Figure [1^ presents the electron density, with two incompressible strips located 
symmetrically near x / d = ±Q.7h and the surrounding compressible regions. With decreasing 
temperature the incompressible strips are stronger pronounced (see inset Fig. [1^) and so are 
the related potential steps (inset Fig.[Tj5). As mentioned in the introduction, in compressible 
regions one of the LLs is pinned at /i*, leading to nearly perfect screening and constant 
potential, while in an incompressible region fi* falls into a gap between the LLs making 
redistribution of electrons, and hence screening, energetically impossible and nei{x) constant. 
We demonstrate these features in Fig.[Tj5, where we exhibit /i* (dash-dotted line), the first two 
LLs (dashed lines with opaque circles), and the total potential (solid line). 

In the right panels we show electron density and potential profiles for various modulation 
amplitudes Vq at a low temperature and a high magnetic field, at which the average filling 
factor in the centre region is around 2.5, i.e., the n = 1 LL is pinned to fi*. With increasing 
Vq the external potential (energy) is raised in the centre and lowered near x/d = ±OA, which 
leads to a lowering of nei{x) near x = and an increase near x/d = ±OA. This redistribution 
of charges screens the potential at the low temperature of Fig. [lj:,d so effectively, that the 
resulting weak modulation of the self-consistent potential for Vq/Eq = 0.01 and 0.02 is not 
seen on the scale of Fig. [TJl. At larger modulation, however, the filling factor z/(0) in the 
centre falls below 2, which means that near x = the self-consistent potential must rise so 
much that the lowest LL n = reaches fi* (as is demonstrated for Vq/Eq = 0.07 in Fig.[TJi). 
Therefore, on both sides of the centre region incompressible strips must develop across which 
the self-consistent potential changes by fkUc. 
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Figure 1. Electron densities [(a),(c)] and electrostatic potentials V{x) [(b), (d)] without 
modulation [(a),(b)] and with external modulation potential Vm{x) — Vq cos{qx), q = 2.57r/ci 
[(c), (d)], where Vq/Eq values are given by the legend. Insets in (a) and (c) demonstrate the 
temperature dependence for kT / Eq = 0.1 x lO^^^ (thick solid Hne), 0.2 x 10^** (dashed line), 
and 0.1 X 10"^ (thin soUd line). 



We see from these results that the screening is highly non-linear, since the change of 
the self-consistent potential produced by the harmonic modulation ^ is step-like rather 
than cosine-like. Moreover, for weak modulation Vq/Eq < 0.02 the resulting variation 
AV = V{0) — V{OAd) of the self-consistent potential is so small that it is not detectable 
on the scale of Fig. [ijl, whereas AV^ for Vq/Eq > 0.04 assumes the constant value hwc- 
Apparently this stepwise increase of AV with Vq is a consequence of the pinning of LLs 
to the electrochemical potential (locally perfect screening), which is stronger pronounced at 
lower temperatures. We now consider this effect in more detail. 

3.1. Tuning the modulation strength at fixed magnetic fields 

In Fig.EK we show the "variance" of the screened potential as function of the amplitude Vq of 
the applied modulation for several values of the magnetic field and a fixed low temperature, 
and in Fig. |2j5 the same dependence for a fixed value of the cyclotron energy and several 
temperatures. We define this variance as the difference between the maximum (at a; = 0) 
and the minimum (near x = 0.4c/) of the difference V{x; Vq) — V{x; 0) of the self-consistent 
potentials calculated with and without modulation, respectively. 

For a fixed value of the magnetic field (i.e., of = hcUc), we obtain a stepwise increasing 
curve, as we expect from the discussion of Fig.[T]:,d, which corresponds to = ^12- The 
corresponding (thick solid) curve in Fig.|2^ shows two broadened steps around Vq/Eq = 0.025 
and 0.075. We have already discussed the reason for the step at the lower Vq value and for 
the plateaus below and above this step. With further increasing Vq the filling factor near the 
density maxima becomes larger than 4 (see Fig. [l]:) while the density in the centre is still 
finite. Then new incompressible strips with local filling u{x) = 4 develop on both sides 
of the density maxima, which corresponds to a decrease of the self-consistent potential by 




Figure 2. (a) Variance of screened potential versus amplitude Vo of the external potential 
for various magnetic field values (Trao/d — 0.01, kBT/Eo — 2 • lO"'"^). The dash-dotted 
line indicates linear screening for zero B and T, see text, (b) Temperature dependence of the 
variance for il/f72 = 1 



the amount to new local minima near x = ±OAd. In Fig. |2k this leads to the plateau 
of height 2^2. As Vq increases further, the electron density nei{0) in the centre vanishes 
and no further redistribution of electrons from the centre region to the region of maximum 
density is possible. Then the high screening ability of the electron system breaks down 
and the variance increases with further increasing Vq much more rapidly (see behaviour at 
Vo/Eq > 0.125). This breakdown happens for Q = ^2 when the filling factor at the density 
maxima is somewhat larger than 4 (~ 5). 

We can now understand all the curves in Fig. |2^. Since, apart from the fine structure 
near incompressible strips, the global density profile depends only weakly on the magnetic 
field, for the lowest curve (long-dashed, = 8^2) typical filling factors i'(x) are by a 
factor 3 smaller than for the case Q = we have just discussed. Therefore the breakdown 
situation is reached when the filling factor near the density maxima is of the order 5/3 < 2. 
Thus the breakdown sets in before an incompressible strip and thus a plateau of height 3^2 
of the variance can develop. For Q = 2^2 [short-dashed line, label (i)] the local filling 
factors z/(x) are typically about half as large as in the case 17 = ^2- Thus, z/(x) < 2 for 
weak modulation, and incompressible strips with ^{x) = 2 can occur only near the density 
maxima for sufficiently strong modulation, and a single plateau of height 21^2 occurs below 
the breakdown regime. A similar situation is met for n = 1.5^2 [label (ii)], but since 
now i'(x) is about 2/3 times that for = 0.2, the threshold of the plateau is reached at a 
considerably weaker modulation. Finally, for f2 = 0.75^22 we obtain three plateaus in the 
pre-breakdown regime, corresponding to the successive development of incompressible strips 
with local filling factor 4 (maxima), 2 (minimum), and 6 (maxima). 

We should mention that the behaviour of the variance curves is not always that regular. 
For special values of it may happen that incompressible strips near the centre and near the 
density maxima occur at the same modulation strength. Then steps of height occur. We 
met such a situation for = 0.51)2 (not shown in Fig.|2^), which yields nearly the same trace 
as r2 = ^l2. 

For the smaller Q values, the variance curves oscillate around the dash-dotted line in 
Fig. EK, which represents the corresponding result 2Vo/[EQ6{q)] of linear Thomas-Fermi 
screening ||T|| for zero T and B, with e{q) = 1 + l/(aog) = 41 for q = 2.5Ti/d. This linear 
screening approximation breaks down when the amplitude of the screened potential (i.e. half 
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the variance) becomes equal to the Fermi energy of the unmodulated system Q, which can 
be estimated from Fig.[l]as ~ 3 x 1{)^'^Eq, corresponding to Vq ~ 0.12_E'o- Finally we should 
mention that the small-variance regime in Fig. |2h also corresponds to linear screening, but 
now with e{q) = 1 + (Dt/Dq) / (aoq), where Dt is the thermodynamic DOS in a LL, which 
can be estimated by Dt/Dq ~ hwd {^kBT) lU]. Thus, the slope of the variance curves is by 
a factor ~ AksT /Vt <^ 1 smaller than that of the dash-dotted line. 

As shown in Fig.|2b, with increasing temperature the plateaus of the variance- versus- Vq 
curves get a finite slope and the steps are smeared out. This is an immediate consequence of 
the temperature dependence of the incompressible strips (see Fig.[Tt). At the highest shown 
temperature screening becomes again linear and independent of B, but now the temperature 
is so high that the thermodynamic DOS is noticeably smaller than Dq. A clear indication of 
steps is seen only for ksT /fvjjc < 1/25 Q. 

3.2. Sweeping the magnetic field at fixed modulation strength 

From the above considerations we expect that, for fixed modulation strength Vq, the variance 
as a function of the cyclotron energy Vt will mainly follow integer multiples of Vt with abrupt 
transitions between neighbouring integers corresponding to the steps in Fig.|2^. These abrupt 
transitions will be similar to the familiar magnetic-field-dependent jumps of the chemical 
potential of a homogeneous 2DEG at fixed density, which have nothing to do with the periodic 
modulation and should be separated from the modulation-induced structures. To this end, we 
first discuss the chemical potential. 

In a Hall bar with spatially varying electron density and electrostatic fields the 
electrochemical potential fi* is a thermodynamic quantity being constant in equilibrium. As 
an analogue of the chemical potential of a homogeneous 2DEG we define = fi* — V{x) 
and denote in the following /i(0) as the chemical potential. In Fig.|3^ we show, together with 
the LLs Vt{n + 1/2), a plot of //(O) versus Vt, which exhibits a saw-tooth shape as known 
from the unbounded 2DEG. To understand this, we plot in Fig.|3j5 and|3J: filling factor and 
self-consistent potential V{x) for four magnetic field values marked in Fig.|3^. For the lowest 
selected B value [opaque circle in Fig.|3K, label (i)], in the centre region the LL with n = 1 
is pinned to /i*, the filling in the centre is z/(0) > 2, and the incompressible strips with 
z/(x) = 2 are close to the edges. As Vt increases these strips move towards the centre, and 
the compressible centre region, with u{x) > 2 and a flat potential minimum, shrinks [case 
(ii), filled square]. In the transition region at still larger [(iii), opaque diamond] the strips 
merge to an incompressible centre region, where the potential minimum is no longer flat, and 
its depth (measured from the adjacent compressible strips) is smaller that fi. As Vl sweeps 
through the transition region, the width of the incompressible centre region and the depth of 
the potential minimum shrink to zero. When the transition is completed [(iv), filled triangle], 
we have again a broad compressible centre region in which the next lower LL (here n = 0) 
is pinned to n*. Thus, in the confined 2DEG in a Hall bar, the jump of the chemical potential 
from a LL to the next lower one is realised by a drastic change of the position dependence of 
the self-consistent potential, which in our idealised model happens in the centre. 

To separate these effects from the modulation-induced screening effects, we define now 
the variance of the self-consistent potential as var = V{0) — V{xrnin), where zLxmin ~ ±0.4(i 
are the positions of the minima of V{x) in the presence of modulation, and V{0) is the local 
maximum at the centre. In Fig.|4^ we plot var versus fl, together with /i(0). As expected, var 
follows widely integer multiples of fi, whereas yu(0) follows half-integer ones. But, whereas 
/i(0) jumps with increasing Q always to the next lower LL, var can also jump to a higher 
multiple of To understand this, we show in Fig.|4j)-|4^ density and potential profiles for the 
six Q values marked in Fig.H^. 
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Figure 3. Chemical potential and LLs (a); selected density (b) and potential (c) profiles. 



In case (i) (solid lines in Fig. EJ? and HJ:) the centre region is compressible with z/(0) 
slightly larger than 2 and with the n = 1 LL pinned to n*, whereas the compressible regions 
near the density maxima have iy{x) > 4 (^{xmin) ~ 5) and the n = 2 LL is pinned to /i*. 
This yields var ^ Vl. For slightly larger [(ii), dashed lines inHJ? andHJ:], an incompressible 
strip with filling factor 2 and a local maximum of V{x)a.tx = develops in the centre, while 
the situation near the density maxima is unchanged. This leads to i7 < var < 2^7. At still 
slightly larger Q [(iii), dotted lines in|3j5 andHJ:], the region near the density maxima is still 
qualitatively unchanged, whereas in the centre a compressible region with z/(x) < 2 develops, 
where the n = LL is pinned to n*, leading to var ^ 2^7. As ^7 increases further [(iv) - 
(vi)], the situation in the centre remains qualitatively the same, whereas the filling factor at 
the density maxima changes from y{xmin) > 4 to ^{xmin) < 4, so that the incompressible 
strips with z/(x) = 4 merge [(v), dashed lines in|4|l and|4^] and finally disappear, accompanied 
by the disappearance of the potential step across these strips. Thus, we have again var ^ fi, 
until for much larger f2 the maximum filling factor becomes smaller than 2 and var goes back 
to the very small value corresponding to linear screening in the lowest Landau level. 

In summary, the low-temperature screening properties of a confined inhomogeneous 
2DEG in strong magnetic fields can be understood from a few basic facts: Incompressible 
regions with integer values of the local filling factor (even integer values for a spin-degenerate 
2DEG), and thus constant electron density, are accompanied by steps of the self-consistent 
potential between the adjacent compressible regions, where adjacent Landau levels are pinned 
to the constant electrochemical potential, so that the step height is JtlUc, the cyclotron energy. 
The modulation-induced variation of the self-consistent potential ("variance") is, therefore, 
usually an integer multiple of hiOc and changes by ±hLUc if incompressible regions appear 
or disappear (near local maxima or minima of the electron density) due to a variation of 
modulation strength or magnetic field. Screening of the applied modulation potential breaks 
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down if the electron density near the density minima becomes so small that no further 
redistribution of electrons from the local potential maxima to the local potential minima is 
possible. 
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